Conductivity in the Heisenberg chain with next to nearest neighbor interaction 
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We consider a spin chain given by the XXZ model with a weak next to nearest neighbor pertur- 
bation which breaks its exact integrability. We prove that such system has an ideal metallic behavior 
(infinite conductivity) , by rigorously establishing strict lower bounds on the zero temperature Drude 
weight which are strictly positive. The proof is based on Exact Renormalization Group methods 
allowing to prove the convergence of the expansions and to fully take into account the irrelevant 
terms, which play an essential role in ensuring the correct lattice symmetries. We also prove that 
the Drude weight verifies the same parameter-free relations as in the absence of the integrability 
breaking perturbation. 

PACS numbers: 05.60Gg, 05.10.Cc, 75.10.Jm 
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INTRODUCTION 

The conductivity properties of quantum spin chains 
has been the subject of an intense research in recent 
times, see e.g. [1-16], using both numerical and analyt- 
ical methods; however several basic aspects remain still 
controversial and this makes the interpretation of exper- 
iments [17, 18] problematic. A prominent role among 
systems modeling spin chains is played by the XXZ 
model, which was solved long ago by Bethe ansatz [19] 
and describes spins with nearest neighbor perturbation. 
If S% = for i = 1, 2, . . . , L and a = 1, 2, 3, cr£ being 
the Pauli matrices, the Hamiltonian of the XXZ chain is 
Ho = 



L-l 



JSlSl +1 + J 3 SlSl +1 + hSl}+UL (1) 



and Ul takes into account boundary conditions. We will 
choose J — 1 for dcfinitcncss. The above model can be 
rewritten as a many body system of interacting spinless 
fermions through the Jordan- Wigner transformation (see 
below); when J3 = the system is equivalent to model 
of free fermions. 

An important question to be understood is how much 
the solvability property of the XXZ chain influences 
the conductivity. Indeed in classical dynamics trans- 
port properties are very much affected by integrability 
or nearly integrability, see e.g. [20, 21], and one could 
wonder if the same happens in the quantum case. We 
can therefore add to the XXZ chain a next to nearest 
neighbor interaction breaking exact solvability, that is we 
consider the Hamiltonian 



H = H Q + H 1 
introduced in [22] , with 

L-l 



Hi — A [S x S x+2 + S x S x+2 + S X S X 



+2J 



(2) 



(3) 



Using the Peierls substitution we can derive the expres- 
sion of the spin current, given by 



jx — S X S X+1 S X S X+1 + XF X 



(4) 



where F x is an expression quartic in the spin operators 
whose explicit form will be written in the following sec- 
tion. If p x = S% + \ and (j x ,jl) = (p x ,j x ) we define 



dx^e' 



x ,pJx ,-p 



>P,T (5) 



and < O > 



Tre 



'O 



Tre-^ , O xo = e Hx «Oe~ Hx <> and T de- 
notes truncation. We will use the notation lim^oo < 
O >p=< O >. An important thermodynamic quantity 
is the susceptibility defined as k\ — lim p _j.o lim po ^o < 
PpP-p > with p = (po,p). 

The framework for most transport studies is linear re- 
sponse theory, where the conductivities are given in terms 
of dynamic correlations. According to Kubo formula, the 
spin conductivity at frequency u> at zero temperature is 
given by 



cr\(u) = lim lira lim — 1 

s^o p->q @->oo ip a 



u+iS 



where p = (po,p) and 



D^ x (p) = [Kl\(p)+ < 3 U > p ] 



(6) 



(7) 



where j D is the diamagnetic term (see below) and D\ = 
lim po _j.o lim p ^o lim^oo Dp^p) is the zero temperature 
Drude weight [23]. At zero temperature a non van- 
ishing Drude weight signals an ideal metallic behavior 
with infinite conductivity; this is what happens at the 
free fermion point J3 = A = 0. In the XXZ chain 
(J3 7^ 0, A = 0), corresponding to interacting fermions, 
the Drude weight can be computed by Bethe ansatz 
[19], [24], [25] and a non vanishing result is also found 



x=l 



7T sm /1 
p 2p(n — p) 



cos p 



■Ja 



(8) 
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implying infinite zero temperature conductivity; this 
should be compared with other interacting ID Fermi sys- 
tems, like the Hubbard model at half filling, in which even 
an arbitrarily weak interaction induces an insulating be- 
havior [26]. Remarkably the XXZ Drude weight is non 
vanishing even at positive temperature, see [1] (h ^ 0) 
and [15] (h = 0); the conserved quantities in the inte- 
grable XXZ chain imply dissipationless current at finite 
temperature. 

Much less is known about the the conductivity in pres- 
ence of an integrability breaking terms as (3). There is 
no consensus even on the basic question if integrability 
breaking terms make the conductivity finite or not at non 
zero temperature; some groups have results supporting a 
finite conductivity [1, 2, 5, 6, 9, 11] while others get ev- 
idence of an infinite conductivity [3, 7, 8, 12, 16]. The 
reason of this ambiguity is the subtle interplay between 
dangerously irrelevant terms and conserved quantities, 
making the results particularly sensitive to regulariza- 
tions or approximations. 

All the technical problems making the understanding 
of the finite temperature conductivity properties of the 
non integrable chain (2) so difficult and the conclusions 
so uncertain appear already at zero temperature; there 
exist indeed no rigorous results on the conductivity for 
the non integrable spin chain (2) even at T = 0. Numer- 
ical analysis using Exact Diagonalization or Montecarlo 
have severe limitations due to the finite system size. On 
the other hand the Bcthe ansatz solution does not furnish 
a good starting point for a perturbative analysis. Usually 
the zero temperature Drude weight is computed replacing 
the chain model by a continuum Quantum Field Theory 
model (the Thirring or the Luttinger model) which can 
be mapped in a boson gaussian model (bosonization) in 
which the Drude weight can be explicitly computed. The 
difference between the continuum and the original lattice 
model is in terms which are irrelevant in the Renormal- 
ization Group sense, which take into account lattice ef- 
fects like Umklapp and non linear bands. Such irrelevant 
terms are crucial ones at T > (neglecting them the 
conductivity is infinite and it has been conjectured that 
their presence can render the conductivity finite) but also 
at T = they are important and cannot be neglected; 
their contribution to the Drude weight is of the same or- 
der than the value found in the gaussian model, so that 
in order to get a lower positive bound for D\ one has to 
exclude cancellations. More in general, in conductivity 
problems one cannot trivially appeal to an universality 
principle and results are strongly sensitive to the irrele- 
vant terms, that is to the details of the model. An exam- 
ple in which this is particularly transparent is given by a 
recent theorem [29] proving the universality of graphene 
conductivity in presence of interaction, in which the role 
if irrelevant terms is crucial. 

Our main result is the following theorem 



Theorem. There exists e < 1 such that, if\J^\,\X\ < e 
the zero temperature Drude weight is non vanishing and 
analytic in J3,A; moreover 



n ts Vs ^ K 
D\ = K^- K\ = 

TT TTV S \ 



(9) 



with K 



1 - 



TTV S \ 



[( J 3 + 2A)(1 - cos 2 PF ) + A(l - cos Ap F ) + F] 



and v s = sin(p F ) + F, s\np F = h and \F\ < Ce 2 , \F\ < 
Ce. 

The above result rigorously establish for the first time 
that the Drude weight is finite for the anisotropic XXZ 
chain perturbed by a weak next to nearest neighbor in- 
teraction (3), so that the system behave at zero tempera- 
ture as an ideal metal (infinite conductivity). It is based 
on the techniques introduced in [31] and it extends [33] 
which was limited to a very special anisotropic spin in- 
teraction. Besides providing a strict lower bound on the 
Drude weight, the above theorem provides the validity of 
the relation (following from (9)) 



Dx 



(10) 



which is known to be true at A = from Bethe ansatz 
[19, 24, 25] as 

= TT1 rrsin^ u s0 = -sin/z (11) 

2tt (it — (i) \i 

Note that for A = h = and small J 3 R- 1 = 2(1 - 
f ) = K- 1 = 1 + 2ii + and Vs = 1 + 0(J 3 ) in 

agreement with our formulas. The validity of (9) at zero 
temperature for the spin chain (2) was conjectured in 
[27], but its validity was only checked at A = when the 
exact solution is valid. 

The proof of the theorem is based on Exact Renor- 
malization Group (ERG) methods, sec e.g. [30], which 
appears quite well suited for the problem; contrary to 
the usual field theoretical RG, in which the irrelevant 
terms are simply neglected, in ERG no approximations 
are done and all the irrelevant terms, which are crucial 
in conductivity properties, are fully taken into account. 



LATTICE WARD IDENTITIES 

The spin chain Hamiltonian can be rewritten as 
a fermionic Hamiltonian by using the Jordan- Wigner 
transformation; calling S*+ = S^ + iS 2 and S~ = S].—iS 2 
we can write 



=- i7r Ei 



(12) 
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and S x = a x a x — 1/2, where a x are fermionic creation 
and annihilation operators. One finds 



L-l 



x = l 



x=l 



-J 3 (ata x - ^)(a+ +1 a x+1 - -) 



and 



L-l 



H i = - x ^[- a t{ a t+i a x+i ~ \) a x+2 



(13) 



(14) 



a x+2( a x + l a x + l ) a x + ( a i 



)(dt+2 a x+2 ~ o)] 



and Ul can be chosen so that a x verify periodic bound- 
ary conditions. The above representation makes straight- 
forward the computation of the spin currents which are 
e.m. currents in the fermionic representation. Using the 
Peierls substitution, the coupling with a classical e.m. 
field is described by the Hamiltonian H + H(A) with 



L-l 



H{A) = -J2 \\ a tUx,x+ia x+1 + a+ +1 U x+hx a x } + 
x=i 1 

L-l 1 

^^2iat (a++i a x+i - ~)Ux,x+2a x+2 + (15) 



x = l 



a x+2( a x+l a x+l - -j) U x+2,xa x ] 

and U Xj y — e' l f* dsA i( s ) — 1. The paramagnetic current is, 
lf Vi,p = — §T~ and m , P = — 



3p = -- 



dH(A) 1 



dA z 



\ -2i{k+p) _ & 2ik 



2l Ath+p + 



-i{k+p) _ e ik 

— 2i — y a iH + p 

(16) 



-Jfi V2,pw(k,p)a+ i a+ 2 a k3 a ki+k2 _ k3+p 

ki-k 2 ,k 3 

With w(k, P ) = ^ e -^-ik 2 +ik 3 -2ip _ e +2ik 1+ ik 2 -ik 3 ^ 

Similarly the diamagnetic current is defined as 



8 2 H(A) 
dApdA^p 1 



(17) 



Note that the current is not anymore quadratic as in the 
A = case. 

The vertex and the two point correlation are connected 
by the following Ward Identities 



-ipo < P P a k d+ +p >f3, T +P < jpK^i+p >P,t= 



[(«k«k> 



fi,T \ a k+p a k+ P/ T . 



(18) 



Similarly the current-current and density-density corre- 
lations obey to 

-^;° a (p)+^ a (p) = o 

-ipotyfa) +P(^;i(p)+ < j D > P ) = o (i9) 

Let us consider the (3 — » oo limit. If Kqo(po, 0) and D(p) 
were continuous in p = 0, (19) would imply that both 
k and D are vanishing. In the case we are considering, 
we will see in the next section that K oo (po,0) and D(p) 
are bounded but not continuous in p = 0, and this fact 
implies the following identities 

X$°(po.O) = 0, [K 1 x 1 {v)+< ] D >] PO=0 = (20) 

The simplest derivation of the above Ward Identities is 
through to the Grassmann integral representation for the 
correlations. We introduce the generating function 



WW = J P(#)< 



,V(V) + (V» + B(A^) 



(21) 



where, if k = (&o,fc) with ko the Matsubara frequency, 
P(dip) is the fermionic gaussian integration with propa- 
gator 



ff(k) = 



— ifc + cos k — [I 
with \i = h + A + J3 and 



v= 

Jo 



x ■f y x+2e 2 



+ i>. 



x+2e 2 



(22) 



C) + (23) 



Js^^x^x+e^x+e! ~ ^^x+ ei ^x+ ei V' x +2e 1 ~ 
Al /'x+2e 1 V4+e 1 V'x+e 1 V'x + ^+Vx V>x+2 ei ^x+2 ei ] 

Finally the source term is given by 

ff> r 

J o x L 

[V ,+ (e '/: +1 ^W - l)C +ei +^ +ei (e 2/ ^ Al(s) - l)Vx] + 

A [ P d X0 J2 X ^(e f:+2Al(S) ~ 1 )V'x+2e 1 

+ ^+ + 2e 1 (^ / ^' 4l(S) -l)^x + (24) 

A(e-'^ +2 ^W - l)^+ +2ei ^ +ei ^- +ei Cl (25) 



and 



/^oB^ + ^xl (26) 
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The correlations are easily written in terms of derivatives 
of the generating function; in particular 

^ d^W{A,<f) 



~ , d*W{A,<p) | 

K ^ (P) = dA , p dA J° 

and so on. Performing the phase transformation 



(27) 



(28) 



in(21), we find 

W(,4 + <9a,0e M ) = W(A,<P) , (29) 

Therefore by performing derivatives with respect to a 
and to the external fields A, 4> the Ward Identities 
(18), (19) follow. 

EXACT RENORMALIZATION GROUP 
ANALYSIS 

The perturbation theory for the correlation functions is 
(apparently) affected by infrared divergences, related to 
the divergence of the free propagator Eq.(22) at cosfc = 
ji. As the interaction modifies in general the location of 
the singularity it is convenient to write u = cos pf + v, 
where v is a counterterm fixed so that the singularity of 
the two point function is at k = pp. We consider the 
following equivalent generating functional 

e w(A,4>) = J p We vW)+W)+fl(^) (30) 
where P(dip) has now propagator 



.9(k) 



— ik a + ^(cos k — cospp) 



(31) 



where v s = sin pp (I +8), vf = sin pp and in V there are 
two new quadratic term proportional to v, 8; moreover 
B(A, tp) is a source term given by B(A, ip) = J dx[^4°p x + 
^4 1 ix] (we just keep the linear part for the computation 
of the current-current correlation). The source term A p 
is assumed with a compact support and we will choose 
v, 8 as function of J3 and A so that the Fermi point of 
the interacting theory is just pp and the velocity v s . 

The functional integral (21) is perfomed in a multi- 
scale fashion using the following two basic properties of 
Grassmann integration. The first is the addition prop- 
erty, which says that 

J P{dil))Fty) = J P(d^) J P(# (2) )F(^ (1) +V (2) ) 

(32) 



where P(#), P(rf0 (1) ), P(<# (2) ) are Grassmann integra- 
tions with propagators g,g^\g^ with g = + g^ 
and F(ip) is an analytic function in i[>. The other is the 
invariance of the exponential, which says that 



/ 



P(dip)e 



v(*p+<p) = e v'(4>) 



with 



00 



(33) 



(34) 



and £ T are the fermionic truncated expectation. Us- 
ing (32) we can decompose ^ as a sum of indepen- 
dent Grassmann fields living on momentum scales 
k— pp| ~ 2 h , pp = (0, ±pf) with h < a scale label and 
cospf = After the integration of the fields with scales 
> h we rewrite Eq.(30) (setting, for simplicity, 0* = 0) 



as e" 



(<h)\ „V (h) (V> ( - h) )+.B (h) (A,^ h) ) 



(35) 



J «=± 

where ip^, u , w = ±1 is the quasi-particle field at the 
Fermi point p F (with quasi-momentum k' relative to the 
Fermi point p F ) and P<h(i>oj) is a fermionic gaussian 
integration with propagator 



,(<h) 



Xfc(k') 



Zh —iko + lov s sin k! + cospf(cos k — 1) 



(36) 

where Xfi(k') is a cut-off function supported in |k'| < 2 h 
and Zh is the effective wave function renormalization. 
The effective potential \>^- h \ip^- h ^) is a sum of inte- 
grals of monomials in ip(- h ^ of order n multiplied by 
kernels wj^Q (xi, . . . ,x n ); similarly the effective source 
is B( h \A,ip(- h> )) is sum of integrals of monomials with 
n ip~ fields and m A-fields multiplied by kernels W^q . 

The scaling dimension is given by, if WZ4 is the number 
of quartic interactions, n 2 ,A and n^A are the number of 
source terms with two and four fermionic lines and n e 
the number of fermionic external lines 



-2(m 4 + n 2 .A + n^ A - 1) + 

ft 

(2m 4 + 2n 4 ,A + n 2 ,A - = 



(37) 



n 2 ,A 



ne 
2 



Therefore, the marginal terms in the RG sense are those 
with n e = 4, tia = 0, n e = 2, ua = 1, and the rele- 
vant are the ones with n e = 2, ua = 0. All the other 
terms are irrelevant, in particular the terms with six 
or more fermionic fields, corresponding to the effective 
multi-particle scattering terms, or the source terms with 
more than two fermionic fields. The integration of the 



■5 



scale h is done, using the addition property (32), writing 



f n Pzu-M^- 1 

J L0 = ± 



(38) 



CV (h) +TZV (h) +CB {h '> +TZB (h) 



with 1Z = 1 — C and CV^ contains the marginal and 
relevant part of the effective interaction 



where 



CV™ = v h 2 h F v h + 8 h F 5 h + \ h F* (39) 

K = / fl rfa; o 51 H V'x^V'x^ 

2" U=± X 

"2 U)=± £C 

f " = y_£ dx ° ^ ( 4 °) 



2 £ 



while in T^V^^ are all the irrelevant terms. In the same 
way 

B^(A^) = j^Z,, h j * A„(p)^(p) (41) 

With J = E w =±l V'ixC.x, Jl = E w =±l W ^ixC,x- 

Starting from (38), we can integrate the field tp^ using 
(32) and (33) and the procedure can be iterated. 

Each single scale propagator g^(x), given by (36) 
with x/i(k') replaced by //j(k'), a smooth function non 
vanishin only in 2 /l ~ 1 < |k'| < 2 h+1 , verifies the following 
bound, for any integer TV 



iy y 11 ~ 1 + [2 h \x\] 



N 



implying 



|fl (h) (x)| ioo <C2 fc \g (h \ X )\ Ll <C2- h 



(42) 



(43) 



The outcome of the above construction is that the kernels 
Wn^m are expressed as series in the running coupling con- 
stants Vk = (j^fejAfe), k > h. The fermionic expectations 
in (34) are expressed in terms of determinants, and ex- 
panding them one obtains a Feynman graph representa- 
tion for the kernels Wn h m- The Feynman graphs are finite 
uniformly in h (this would be not true in the non renor- 
malized expansion C — 0); however their number grows 
as 0(l\ 2 ), if I is the order, so that in this way can only 
prove that the 1-th order is bounded by C l l\[maxk>h\vk\] 1 
(a result usually called nl bound ), a result which is not 
sufficient to establish the convergence of the series. In 



order to improve such bound one has to notice that the 
fermionic expectation are given by 

£ h (^ 1 ...^+ 1 ..^+J - det^Hx, - y,)) (44) 

and the fermionic propagator can be written in the form 

/■! 

S (/l) (x-y)= / dzoY,AUx-*)-B h {y-z), (45) 



with 



\\M 



dz \A h (z)\ 2 < C2 



-2h 



\B h \\ 2 < C2 



Ah 



(46) 

for a suitable constant C. According to the Gram in- 
equality, if a n x n matrix H has the form Hij = (fj, gj) 
then 



|Det£T| < ||f|| n ||g| 



(47) 



Therefore, while expanding the determinant in (44) and 
bounding each term a bound C n n\2 nh is found, using 
the Gram inequality one gets an estimate without facto- 
rials, namely C n 2 nh . In fact in the kernels not the simple 
expectations £ but the truncated ones £ T appear; how- 
ever one can use the Battle-Brydges-Federbush formula 
allowing to write the truncated expectations as sum over 
chains of propagators (ensuring the connection) times de- 
terminants which can be bounded by the Gram inequal- 
ity. Therefore, see Theorem 3.12 of [31] for the proof, by 
using the above ideas it is proved that the kernels W^L 
are analytic in the running coupling constants (with a 
small but uniform in h radius of convergence) and satisfy 
the bounds 

f \W^%\ < c n+m \\ h \ m 2^- n ' 2 - m ^ h , (48) 

Of course all the consistency of the method rely on the 
fact that the running coupling constants remain in the 
analyticity radius. It has indeed proved in [32] that, by 
suitably choosing v, 6, see [32] 



Ah -> A_ co (A, J 3 ) 5 h ->0,u h ->Q 



(49) 



where A_co(A) arc analytic functions of the coupling; 
in other words, there is a line of fiixed points. This is 
consequence of a property, called vanishing of the Beta 
function, which has been proved in [32] by implementing 
Ward Identities at each Renormalization Group iteration. 

Another crucial property is that the wave function 
renormalization Zh appearing in (38) and the vertex 
renormalizations Zo,h,Zi.h appearing in (41) have an 
anomalous behavior consisting in a power law divergence 
Z h = 0(2^ h ) and Z i>h = 0{2^ h ), with 77,77; positive and 
O(Aq); remarkably such exponents are equal 77 = 7/0 = 771, 
that is 



Z h 



1 + A(A,J 3 ) 



z h 



1 + B(X,J 3 ) (50) 
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and A,B = O(X ), and 77 = aXf t + O(Ajj). The identity of 
such exponents is due to emerging relativistic symmetries 
(see below) which are broken by the irrelevant terms; this 
is reflected by the fact that A ^ B as one can verify 
by an explicit computation. Such emerging relativistic 
symmetries emerges from the following decomposition of 
the fermionic propagator 



no justification in neglecting the irrelevant terms, as they 
give an O(l) contribution to the Drude weight. 
On the other hand xj^ n (x) = 



h=— 00 u>=± 



(51) +cos(2 PF x)[^] 2 5 i';l ? (x) 5 i' l 2^(-x)(l + r, l (x))} 



with 



and 



Ah) 



A(k) 



BL ^ —iko + u;fct; s 

k 



r£ h >(x)|<C 



)2h 



l + [2 fc |x| 



JV 



(52) 



(53) 



for any TV, that is with an extra 2 h with respect to the 
bound for the dominant part. The above decomposition 
says that the single scale propagator is closer and closer 
to the one of a massless Dirac fermions, that is <j^f]j(x), 

plus a correction r£'^(x) taking into account of the non 
linear bands. Of course even if the relative size of the 
two terms in (51) is asymptotically vanishing, the con- 
tribution of the r£f^(x) terms to the thermodynamical 
constants is not negligible. The current current correla- 
tion is naturally decomposed in two terms; one, denoted 
by K^ 1X {-x), which takes contributions from the non ir- 
relevant part of the effective potential CV^ and from 
the "relativistic" part of the propagator g^ R (x) plus a 
rest -ftr^ n (x) depending from the irrelevant terms 



(°)H/ \ 1 7^(6)11/ \ 
(x)+K\' (x) 



where 



*? )U (x) 



E 

h— — 00 



s£>(x) 



and 



\sf{ 



[2"|x|] 



N 



(54) 



(55) 



(56) 



where the extra 2 , 6 = 1/2, with respect to the dimen- 
sional bound is due to the fact that it has contribution 
from the irrelevant terms. From (56) and (55) we obtain 
the following bound 



\K 



(6)11 



(x)|< 



C 



1+ x 



2+6 



(57) 



from which we deduce that the Fourier transform 
K^ n (p) is continuous in p and 0(1). There is therefore 



where |0T fc (x)|, |0"f h (x))| < C\X \2- nh , and Z %h ~ 
2^ with fj = bX + O(Xl). We are interested to the 
Fourier transform K^ 11 (p); the contribution from the 
oscillating part of (58) is surely bounded close to p = 
but it is not obvious at all that Fourier transform of the 
first term (the non oscillating one) is bounded. It be- 
haves for large distances, from (50), as 0(|x|~ 2 ) and log- 
arithmic divergences could be present. In conclusion, the 
exact Rcnormalization Group analysis provides non per- 
turbative bounds for the current-current correlations in 
the coordinate space; we cannot deduce by such bounds 
even that the Drude weight is finite, as logarithmic di- 
vergences could be present. Even if we could prove that 
K^ n (p) is bounded in the limit this would not allow us 
to conclude anything on the Drude weight for the pres- 
ence of £f )n (p) which is 0(1). 



EMERGING CHIRAL SYMMETRIES 

The dimensional bounds obtained from the Renormal- 
ization Group analysis are not sufficient for the Drude 
weight and one needs to exploit both the lattice and the 
emerging chiral symmetries of the theory. In order to do 
that, and remembering (51), we introduce a model ex- 
pressed of massless Dirac fermions in d = 1 + 1 with light 
velocity v s . The generating functional is given by 



W rel (B,, 



-I 



P(#)e 



V(v / Z^) + (i/-,0)+B(B,V) 



(59) 



where P{dip) is the fermionic integration with propagator 

1 XJv(k) 



5a; (k) 



Z —ik + Lov s k 



(60) 



with Xiv(k) a cut-off function selecting moments less than 

2 N 



V = X oa J dx J dyv(x - y)^X^+- 
1 

B( J B,V) = E B ^ Wj ^ 



(61) 



with v(x) a short ranged interaction \v(p)\ < Ce K ' p ', 
^ = «+^x,+ + <-V>x,-) and Ji = «+Vx,+ - 
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?/>x _?/> x - )■ Note that in this case x is a continuum vari- 
able both in space and time while in the previous case 
the the spatial component was discrete. The Schwinger 
functions are given by derivatives of the generating func- 
tion 



G^(k,k + p) 



d 3 W, 



rel 



i = 1,2 



d 2 W rel 



AO 2 l \ u vv re i 



(62) 



We can analyze the above functional integral using a 
Renormalization Group analysis. With respect to the 
previous case, in which the lattice furnishes an ultravi- 
olet cut-off and the problem is an infrared one (that is 
the zero temperature and infinite volume limit), in the 
model (59) there is both an ultraviolet and infrared prob- 
lem. We write then -0 = X)-oo where ip^ lives on 
momentum scale |k| ~ 2 h . Note the crucial difference 
with respect to the chain model described above; in the 
present case there are positive scales (the momentum k 
is unbounded in the N — > oo) while in the presence of 
the lattice the scale are h < 1 as k e [— n, it] for the pres- 
ence of the lattice (there is only an innocuous ultravio- 
let problem for the unboundedness of ho). The integra- 
tion of the positive ultraviolet scales has been analyzed 
in [34], and it has proved that after the integration of 
the fields ^( JV ),^;( JV_1 ), one gets an effective po- 

tential with kernels uniformly bounded as N — > oo. A 
crucial role in establishing this result relies on the non- 
locality of the interaction, which eliminates the possible 
ultraviolet divergences present in the case of local 5- like 
interactions. Once that the ultraviolet scales are inte- 
grated out, the integration of the negative infrared scales 
is done as described above for the spin chain, with some 
obvious modification due to symmetry; we call the cor- 
responding effective couplings Xh, Zq^-, Z\ t h while due to 
symmetry 5^ = 1^ = 0. Again the beta function is van- 
ishing and the effective coupling tends to a line of fixed 
points Xh — >h^-oo A_oo, with A_oo an analytic function 

of Aoo. 

The model (59) can be considered the continuum limit 
of the chain model. It verifies more symmetries than the 
chain model; for instance it is symmetric with respect 
to space-time inversion and it invariant under the chiral 
transformation ^ x — > e ±,au ^ x . The advantage of our 
ERG method is that the relation between the two models 
can be understood quantitatively in the following precise 
sense: it is possible to choose Aoo, Z, Z^°\ Z^ so that, 
for small p, k', lo = ± 



with i?i,i? 2 = 0(|k'|,p) and 

K^(p) = G°£(p)+A^(p) 



(64) 



with Ap iV (x) < C|x|~ 2 ~ 9 , hence A M)I/ (p) is contin- 
uous in p at p = 0. This can be proved choos- 
ing Aqo, Z, Z(°\ Z^\ by the implicit function theorem, 
so that the differences between \h, Zh, Zo,h, Zi,h an d 
Xh, Zh, Zo,h, Z\_h is asymptotically vanishing as h — > — oo 
as 0(s2 8h ) from (51). It turns out that 

Aoo = 2(J 3 + 2A)(l-cos2p F ) + 2A(l-cos4p F ) + f (65) 

with \F\ < Ce 2 , and that Z ( 0) ^ Z^ (the symmetry 
between space and time is broken in the chain model). 
In order to understand how (64) is derived we can sim- 
ply notice that if^ n (x) is identical to G"^(p) replac- 
ing Xj, Z h , Z a ^ h , Zi t h with Xh, Z h , Z ,h, Z\ t h\ to the func- 
tions A li:U (p) contribute terms coming from such differ- 
ence and from A^ n (x) in (55). They are all bounded 
by (56) summed over the scale and therefore the Fourier 
transform is continuous in p. The r.h.s. of (64) says 
that the current-current correlation of the chain model is 
equal to the one of a continuous relativistic model, but 
the contribution of the irrelevant terms is not small, but 
simply more regular in Fourier space. 

As we noticed in the previous section, from the dimen- 
sional bound we cannot exclude that the Fourier trans- 
form current-current correlation of the chain model has 
logarithmic divergences; this can be achieved by using 
(64) which allow us to exploit the symmetries of the 
model (59). Indeed while in the spin chain model there is 
only one set of Ward Identities, in the model (59) there 
are two set of ward identities, related to the global and 
chiral symmetry; by performing the phase transforma- 
tions V'w.x — > e mx (global phase transformation) one 
gets, as proved in [33] 

zHp * k + p) + p« * £*;i(k, k + P )] = 



z(°) ZQ 

= A[G£°(k)-G£°(k + p)] 



while performing the phase transformations 
e lu,a *Vw,x (chiral phase transformation) it is found 



(66) 



Z[-ip J^G^(K k + p) + v s p ^yGo;i(k, k + p)] 
= ojA[G 2 f(k)-G 2 f(k + p)} , 
with 



< a k'+p£ a k'+ P 



>=<#°(k')(l + 0(k')) (63) 

< P P «k' +P £«k< +P+P £ >t= Goi(k',k' + p)(l + iZx) 

< 3p«k' +P ^k' +P+P ^ >t= G 2 £0d, k' + p))(l + R 2 ) 



1 - r, A~ 



1+T 



4ttv s 



(67) 



Note the presence of r in the above Ward-Identities, 
which represent the chiral anomaly. In the same way 
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there are two set of Ward Identities for the densities, re- 
lated to the global and chiral transformations; from them 
one can write an explicit expressions the density correla- 
tions, that is, if D u (p) = — ipo + uv s p 



G 



1 (Z(°)) 2 D ( p ) , D+(p) 



0,0 



4tto s Z 2 1-t 2 L L>+(p) D-(p) 



2t1 



1 (Z^) 2 £>_(p) fl+(p) _ 
4™ S Z 2 1 - r 2 [ £ + (p) + £)_(p) Zn ™> 

From the above expression is easy to verify that G°' 2 (p) 
are no£ continuous in p. The density and current cor- 
relations are symmetric between exchange of space and 
time, contrary to what is expected for the chain model 
in which the space and time symmetry is broken by the 
irrelevant terms. Note also that Ward Identities alone 
allow to get an explicit form for the density correlations 
in the effective model (59), what is not possible in the 
chain model (2). 

Eq. (64) provide a relation between the current-current 
correlation of the chain model (2) and of the rela- 
tivistic model (59), once that its bare parameters Aoo, 
are properly fine tuned. The WI (66) com- 
bined with (64) must coincide with the lattice Ward Iden- 
tities (18); this provides relations between the bare pa- 
rameters, namely 

^ - 1 (69) 
Z« 



1 



z(°) _ 1 



Let us consider now (64) in which by continuity 
linip^o ^^(p) = -4pu/(0) (the limit does not depend from 
the order contrary to what happens in the first term in 
the r.h.s. of (64)). The value of the constant A^ V {Q) is a 
complicate function depending from all the micro detail 
of the chain model; however its value is fixed by (20) 



Goo 2 bo,0) + ioobo,0) =0 
G°^(0,p) + A u (0,p)+ <j D >=0 



(70) 



so that 



A) (0,0) : 

A n (0,0)+<j D > 



lim Urn G° g(p) 

p ->-0p->0 



■ lim lim G?f(p) (71) 

p— >0 po-+0 



Note that Aqo(0,0) and An(0,0), containing the contri- 
butions from the irrelevant terms, are essential for the 
result; they break the symmetry between space and time 
which is present in the effective model (59) which is not 
true in the chain model (2); as a result the limits limp-^o 
and lim po ^j.o do not commute in the current correlations. 



Therefore by (64), (68) ,(71) 



DxiP) = 



1 



(Z( 



0)\2 



2 2 
V sP 



irv s Z 2 1 
1 (Z«) 2 



r 2 p 2 



v 2p2 



vl 



TTVsZ 2 1 



T 2 pi 



v 2p2 



O(p) (72) 



O(p) (73) 



and using (69) 

*T(p) - 

Dx(p) = 



K v 2 p 2 
nv s p 2 + v 2 p 2 
Kv s pi 

7T pl+V 2 p 2 



O(p) 



+ 0( P ) 



(74) 



(75) 



with K = y+7j an d using that K\ = 

lim p ^o lim P0 ^ ^a°(p) and D >- = lim P o^o lim p ^ D x (p) 
finally (9) follows. 



CONCLUSIONS 

We Uave rigorously established a strict lower bound for 
the zero temperature Drude weight which is strictly pos- 
itive for the the XXZ chain perturbed by a weak next 
to nearest neighbor interaction (2), proving that the sys- 
tem has an ideal metallic behavior (infinite conductiv- 
ity). We have also proved that the Drude weight verifies 
the same exact relation (10) as in absence of the inte- 
grability breaking perturbation. The main difficulty in 
the analysis relies in the irrelevant terms, whose role is 
essential in ensuring the correct lattice symmetries, and 
this makes the use of the ERG quite well suited for the 
problem. The results are obtained exploiting both the 
lattice and the emerging chiral symmetries of the theory, 
and relying on rigorous estimates on the large distance 
decay of the correlations, based on determinant bounds. 
Our results are conclusive regarding the properties of the 
Drude weight in a non integrable spin chain at zero tem- 
perature, at least for weak perturbations, and the main 
problem remains the conductivity properties at non zero 
temperature. The main difficulty of such a problem relies 
on the subtle interplay between dangerously irrelevant 
terms and conserved quantities; therefore we believe that 
the understanding of such an issue at zero temperature, 
which is achieved in the present paper, is an essential pre- 
requisite for an analytical understanding for the positive 
temperature problem. 
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